Analytical Study of Electronic Structure in Armchair Graphene Nanoribbons 
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We present the analytical solution of the wavefunction and energy dispersion of armchair graphene 
nanoribbons (GNRs) based on the tight-binding approximation. By imposing hard-wall boundary 
condition, we find that the wavevector in the confined direction is discretized. This discrete wavevec- 
tor serves as the index of different subbands. Our analytical solutions of wavefunction and associated 
energy dispersion reproduce the numerical tight-binding results and the solutions based on the k ■ p 
approximation. In addition, we also find that all armchair GNRs with edge deformation have energy 
gaps, which agrees with recently reported first-principles calculations. 
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I. INTRODUCTION 

Graphene, as a promising candidate of future nano- 
electronic components, has recently attracted intensive 
research attention ±&2d£JL Graphene consists of a single 
atomic layer of graphite, which can also be viewed as a 
sheet of unrolled carbon nanotube. Several anomalous 
phenomena ranging from half-integer quantum Hall ef- 
fect, non-zero Berry's phased, to minimum conductivity 2 -, 
have been observed in experiments. These unusual 
transport properties may lead to novel applications in 
carbon-based nanoelectronics. In addition, the carri- 
ers in graphene behave as massless relativistic parti- 
cles with an effective 'speed of light' c* ~ 10 6 m/s 
within the low-energy range (e < 0.5eV)£ These mass- 
less Dirac fermions in graphene manifest various quan- 
tum electrodynamics (QED) phenomena in the low en- 
ergy range such as Klein paradox phenomenon. 6 Rib- 
bons with a finite width of graphene, referred to as 
graphene nanoribbons (GNRs), have also been studied 
extensively.^^ ' ^i 11 ! 12 ? 13 ! 14 Recent experiments by us- 
ing the mechanical method 2,3 and the epitaxial growth 
method 4 -*^ show it is possible to make GNRs with vari- 
ous widths. 

The carbon atoms on the edge of GNRs have two 
typical topological shapes, namely armchair and zigzag. 
The analytical wavefunction and energy dispersion of 
zigzag nanoribbons have been derived by several research 
groups * 16 ' 17 For armchair GNRs, the analytical forms of 
wavefunctions within the low energy range have been 
worked out based on the effective-mass approximation. 14 
It is predicted that all zigzag GNRs are metallic with lo- 
calized states on the edges ; 8 ' 9 ' 16 ' 17 while armchair GNRs 
are either metallic or insulating, depending on their 
widths jli^ i 10 ! 14 ! 16 To-date, there is no general expression 
of the wavefunction in armchair GNRs. In this paper, 
we derive a general analytical expression of wavefunction 
and eigen-energy in armchair GNRs applicable to various 
energy ranges. In part II, we focus on perfect armchair 
GNRs without any edge deformation and derive the en- 
ergy dispersion by imposing the hard- wall boundary con- 
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FIG. 1: (color online) Structure of an armchair graphene 
nanoribbon, consisting of sublattices, A and B. The width of 
the armchair GNR is n. Every unit cell contains n number of 
A and B sublattices. Two additional hard walls (j — 0, n+ 1 ) 
are imposed on both edges. 



dition. Due to the quantum confinement, the spectrum 
breaks into a set of subbands and the wavevector along 
the confined direction becomes discretized. We observe 
that the electronic structure of perfect armchair GNRs 
strongly depends on the width of the ribbon. The system, 
for instance, is metallic when n = 3m + 2 and is insulat- 
ing otherwise, where m is an integer ,LSi & 10 i 14 i 16 Further- 
more, we study the low energy electronic structure. The 
linear dispersion relation is observed in armchair GNRs. 
In part III, we evaluate the effect of deformations on the 
edges on the electronic structure of armchair GNRs. Cal- 
culation results based on the derived analytical wavefunc- 
tion show that all armchair GNRs have nonzero energy 
gaps due to the variation of hopping integral near the 
edges. This observation is in line with the recently re- 
ported first-principle calculations^ 
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II. PERFECT ARMCHAIR GRAPHENE 
NANORIBBON 

The structure of armchair GNRs consists two types of 
sublattices A and B as illustrated in Fig. 1. The unit cell 
contains n A-type atoms and n B-type atoms. Based on 
the translational invariance, we choose plane wave basis 
along the x direction. Within the tight-binding model, 
the wavefunctions of A and B sublattices can be written 
as 

1 n 



x Ai 1=1 



' B N B 



(i) 



where 4>A(i) and (j>B (i) are the components for A and B 
sublattices in the y direction, which is perpendicular to 
the edge. and are the wave functions of the p z 
orbit of a carbon atom located at A and B sublattices, 
respectively. To solve </>a(*) and ^b(z), we employ the 
hard-wall boundary condition 



MO) = <M0) = o 

(f> A (n + l) = <t> B {n + l) = 0. 
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(c) n=8 
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Choosing <^(i) = = sin( — and substituting 

them into Eq. (2), we get 
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(3) 



q y is the discretized wavevector in the y direction and 
a = 1.42A is the bond length between carbon atoms. 
To obtain the normalized coefficients, N A and N B , we 
introduce the normalization condition 

a(iI>\ ^)a=b (fl>\fl>) B = l. 



AT x (n+l) 



It is straightforward to obtain Na = Nb = 
where N x is the number of unit cells along the x direction. 
The total wavefunction of the system can be constructed 
by the linear combination of i\>a and t/'b 
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FIG. 2: (color online) Electronic structures of perfect arm- 
chair GNRs with various widths, (a) n=6 (b) 7 and (c) 8, 
respectively. The wavevector is normalized based on the prim- 
itive translation vector of individual GNRs. The value of p 
for each subband is labeled in the figure. 

where denotes the nearest neighbors. 

In perfect armchair GNRs, we set fcjj = t and Si = e. 
By Substituting Eq. (4) and Eq. (5) into the Schrodinger 
equation, we can easily obtain the following matrix ex- 
pression: 



(4) 



e fi 
[i* e 



C A 
C b 



E 



C A 



(6) 



cos(^f^) + e 



-ik x a\ 



where fj, =a (ip\ H \tp) B = t(2e l 

Solving Eq. (6), we get the energy dispersion and wave- 
function as 

E = e±\fi\, 



)± = ^(Ma± ] [^Mb)- 



(7) 



Here, ± denotes the conduction and valance bands, re- 



Under the tight binding approximation, the Hamiltonian spectively. -f < < § is required within the first 

Brillion zone. These results are valid for various energy 
ranges. 

(5) Fig.[2]shows the energy dispersion for perfect armchair 

GNRs with width n = 6, 7 and 8. Here, we set s = 0. 
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FIG. 3: (color online) The first conductance and valence 
bands within the first Brillion zone: exact solutions from Eq. 
(7) (red solid line) and low-energy approximation from Eq. 
(10) (blue dash line) for armchair GNRs with various widths, 
(a) n = 21 (b) n — 22 (c) n = 23, respectively. The wavevec- 
tor is normalized based on the primitive translation vector of 
individual GNRs. 

The results are the same as those obtained by using 
the numerical tight-binding method. The electronic 
structures of armchair GNRs depend strongly on their 
widths. When n — 8, the lowest conduction band and 
the upmost valence band touch at the Dirac point, 
which leads to the metallic behavior of n — 8 armchair 
GNRs. Armchair GNRs, however, are insulating when 
n = 6 and n — 7. Armchair GNRs with the width of 
n = 3m + 2 (m is an integer) are generally metallic 
and otherwise are insulating^* 1 ^ In addition, we observe 
several interesting features in the band structures of 
armchair GNRs. 

(i) If n=7, a flat conduction/ valence band (p = 4) 
exists as shown in Fig. [5] (b). Such a flat band generally 
corresponds to ^-j- = | or equivalently cos = 0. 
The energy dispersion becomes independent of k x and 
the eigen-energy always equals ±|t|. Flat band, in 
general, exists only when n is odd. 

(ii) The subbands can be labeled by the quantum num- 
ber p. Combined with the wave number k x along the x 
direction, the quantum number p can be used to define 
the chirality of the electrons in quasi-lD graphene rib- 
bons similar to that in 2D graphene. To identify differ- 
ent subbands, we need the the quantum number p t of the 
ith conduction/ valence band. Here, the definition of the 



sequence of subbands is referred to as the value of eigen- 
energy, Eq, in the center of first Brillion zone [k x = 0), 



Er 
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2 cos 
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(8) 



For the metallic armchair GNRs with width n — 3m+2 
when — or equivalently p = 2m + 2, the energy 
gap between conduction and valence bands is zero. 
Therefore, p\ = 2m + 2 corresponds to the first conduc- 
tion/valence band in n — 3m + 2 GNRs. For the second 
conduction/ valence bands, Ec should have the minimal 
nonzero value compared to the third or even higher 
band. After analyzing the value of Eq, we find that 
P2 = 2m + 3, P3 — 2m + 1 for metallic armchair GNRs 
(n = 3m + 2). By similar analysis, for n > 10, we can 
obtain p\ = 2m + 1, p% = 2m, p% = 2m + 2 for n = 3m 
armchair GNRs and p\ = 2m + 1, p2 = 2m + 2, = 2m 
for n = 3m + 1 armchair GNRs, respectively. For all 
subbands, there is no general rule of the subband index p. 

(iii) Lots of research interest have been focusing on the 
energy dispersion and wavefunction of 2D graphene and 
ID GNRs within the low-energy range d 11 ' 14 Low-energy 
electrons behaves as massless relativistic particles in a 
2D infinite graphene systemi^^ 1 ^ Whether electrons 
keep their relativistic property when they are confined in 
quasi-lD graphene nanoribbons is an interesting issue. 
In what follows, we will focus on the expansion of our 
analytical expressions to the low energy limit. When 
— > |7r and — > 0, we rewrite the eigenenergy in 
Eq. (7) as 



E « ±§at 



7,2 

A. , . 



±hvpk, 



(9) 



where q y (p) = ~ f 71 "); P ^ s the subband index. 

This low energy expansion generates the E oc k linear 
dispersion, with Fermi velocity vp « 10 6 m/s. This ex- 
pression reproduces the result of k • p approximation. 14 
Note that the wavevector in the confined direction (q y ) is 
discretized, corresponding to different subbands. What is 
worthy mentioning is that this energy dispersion works 
well only at the low-energy limit. By substituting the 
value of pi into Eq. (9), we get the low energy expansion 
of the first conduction/ valence band for armchair GNRs 
as 



£1(3777)^° ± lot J hi + ( 



2n 



3%/3(3777 + l)o 



E 1 (3m+1) K ^° ± \atJkl + { 



^ fex-»o 3at , 

£ x (3m + 2) » ± — k x . 



2tt 



3V5{3m + 2)a' 



(10) 



Fig. 3 shows the quality of low energy approximation. 
For large width armchair GNRs, low energy approxima- 
tion seems work well except the edge of first Brillion Zone. 
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FIG. 4: Local density of the states in the first conduc- 
tion/valence band at k x = for armchair GNRs with var- 
ious widths, (a) n=50 (b) n = 51 and (c) n = 52, respec- 
tively. (These values of n are so chosen to match the results 
in reference^) 



III. ENERGY GAP AND WAVEFUNCTION 
FOR EDGE-DEFORMED GNRS 

Because every atom on the edge has one dangling 
bond unsaturated, the characteristics of the C — C bonds 
at the edges can change GNRs' electronic structure 
dramatically ! 10 i 19 To determine the band gaps of GNRs 
on the scale of nanometer, edge effects should be consid- 
ered carefully. The change of edge bonds length and an- 
gle can lead to considerable variations of electronic struc- 
ture, especially within the low energy range d 1 ' 12 In pre- 
viously reported work, the edge carbon atoms of GNRs 
are all passivated by hydrogen atoms or other kinds of 
atoms/molecuresi 10 ' 11 ! 12 ' 14 ! 19 The bonds between hydro- 
gen and carbon are different from those C — C bonds. 
Accordingly, the transfer integer of the C — H bonds and 
on-site energy of carbon atoms at the edges are expected 
to differ from those in the middle of GNRs. The bond 
lengths between carbon atoms at the edges are predicted 
to vary about 3% — 4% when hydrogeneratedi^i Corre- 
spondingly, the hopping integral increases about 12% ex- 
tracted from the analytical tight binding expression i n i 18 
To evaluate the effect of various kinds of edge deforma- 
tion, we carried out general theoretical calculation and 
analysis with our analytical solution of armchair GNRs. 
In general, we can set the variation of the transfer inte- 
ger and on-site energy as dtij, e, for the ith A-type or 
B-type carbon atom in the unit cell. The Hamiltonian of 
the GNRs with deformation on the edge can be rewritten 
as 



As the width gets larger, the quantum confinement due 
to the edge becomes less important and the ID nanorib- 
bons tend to behave like 2D graphene. For large n, as 
expected, the band structure generates the linear disper- 
sion relationship, E oc k x , in the low-energy limit. 



H = e i '£\i){i\- ]T (t + 5t itj )\i){j\ 



(11) 



<hj> 



It is readily to obtain the energy dispersion and wave- 
function by solving the Schrodinger equation 

£ = 7± \fj, + Sfi\ 



In addition, from the expression of wavefunction, we 
also obtain the local density of electronic states in per- 
fect armchair GNRs, P A (i) = P B (i) oc sin^^i) . Fig. 
4 shows the squared wave functions of the lowest con- 
duction band at the center of first Brillion Zone. Note 
that Fig. 4(a) and (c) reproduce the results of the k • p 
approximation.— The state density oscillates as a func- 
tion of the lattice position. The oscillation period is re- 
lated to S-ti. For n = 3m +2 armchair GNRs, the oscilla- 
p _ 

tion period is just 3, which is shown clearly in Fig. 4(a). 
For n — 3m, 3m+l armchair GNRs, we should write 
into an irreducible form ^. The oscillation period then 
equals a, which is the numerator of the irreducible form 
of To match the results presented ir>ii, we choose 

ii = 51 and n — 52 as an example. We get i — 51 and 
i = 52, respectively. As shown in Fig. 4(b) and (c), the 
oscillation period of state density for n = 51 and n = 52 
armchair GNRs equals their width. 



(12) 



where 7 = ^-j- ^2 £ i sm2 (^pr*) is the energy shift Orig- 
in 1 

inating from the variation of on-site energy, while the 
energy shift from the hopping integral variation is 



n + 



i=l 



l n + 1 



pir 



pre 



sm( — — 1) sin — —(i - l))e 2 
n + 1 n + 1 

+o*t(A)»+i(B) sin(— —i)sm(— —{1 + l))e 2 J. 



"n + 1 n + 1 



(13) 
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Such a general expression could include various kinds 
of edge deformations, ranging from the quantum con- 
finement effect due to the finite width, to the effect of 
saturated atoms or molecules attached to edge carbon 
atoms. This result shows that the deformation leads to 
a considerable deviation of the energy dispersion relation 
and wavefunction of the deformed system from those in 
perfect armchair GNRs. The local state density on both 
kinds of sublattices, however, remains the same as that 
in perfect armchair GNRs. The reason is that the wave- 
functions of sublattices A and B change their relative 
phases, but keep the magnitudes unchanged. The varia- 
tions from both the on-site energy and hopping integral 
contribute to the energy shift, while, the change of on- 
site energy has no contribution to the wavefunction as 
shown in Eq. (8). To verify our findings, we reproduce 
the energy gap observed in the recent workii by consid- 
ering only the variation of hopping integrals of the bonds 
on the edges (Stn = St nn = St, others equals zero). The 
corresponding energy gaps for different width ribbons are 
as follows: 



A3 m +2 



A 



86t 



■ sin 



3m + 1 
A° 

3m+1 3m + 2 



sm 



3m + 1 

2 (TO+ 1)7T 

3m + 2 ' 



= A°. 



2St 



3m+2 m+1' 



(14) 



where A 3m , A 3m+ i and A 3m+2 are the energy gaps of 
perfect armchair GNRs. Their values can be extracted 



from Eq.(8): 2t 



2 cos 



(2ro+l)7T 
3m+l 



+ 1 



2t 



2 cos 



(2m+l)-rr 
3m+2 



1 



and 0. This result suggests all armchair graphene ribbons 
with edge deformation have nonzero energy gaps and are 
insulating correspondingly. 



IV. CONCLUSION 

In this paper, we study the electronic states of arm- 
chair GNRs analytically. By imposing the hard-wall 
boundary condition, we find the analytical solution of 
wavefunction and energy dispersion in armchair GNRs 
based on the tight-binding approximation. Our results 
reproduce the numerical tight-binding calculation results 
and the solutions using the effective-mass approximation. 
We also derive the low-energy approximation of the en- 
ergy dispersion, which matches the exact solution except 
for the edge of first Brillion zone. The linear energy dis- 
persion is observed in armchair GNRs in the low energy 
limit. In addition, we also evaluate the impact of the edge 
deformation on GNRs and derive a general expression of 
wavefunction and energy dispersion. We can reproduce 
the energy gap for hydrogenerated armchair GNRs pre- 
sented ir*ii. When we consider the edge deformation, 
all armchair GNRs have nonzero energy gaps and thus 
are insulting. Overall, the derived analytical form of the 
wavefunction can be used to quantitatively investigate 
and predict various properties in armchair graphene rib- 
bons. 
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